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In  this  work,  a  model-based  fault  diagnosis  methodology  for  PEM  fuel  cell  systems  is  presented.  The 
methodology  is  based  on  computing  residuals,  indicators  that  are  obtained  comparing  measured  inputs 
and  outputs  with  analytical  relationships,  which  are  obtained  by  system  modelling.  The  innovation  of 
this  methodology  is  based  on  the  characterization  of  the  relative  residual  fault  sensitivity.  To  illustrate  the 
results,  a  non-linear  fuel  cell  simulator  proposed  in  the  literature  is  used,  with  modifications,  to  include 
a  set  of  fault  scenarios  proposed  in  this  work.  Finally,  it  is  presented  the  diagnosis  results  corresponding 
to  these  fault  scenarios.  It  is  remarkable  that  with  this  methodology  it  is  possible  to  diagnose  and  isolate 
all  the  faults  in  the  proposed  set  in  contrast  with  other  well  known  methodologies  which  use  the  binary 
signature  matrix  of  analytical  residuals  and  faults. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  energy  generation  systems  based  on  fuel  cells  are  complex 
since  they  involve  thermal,  fluidic  and  electrochemical  phenomena. 
Moreover,  they  need  a  set  of  auxiliary  elements  (valves,  compres¬ 
sor,  sensors,  regulators,  etc.)  to  make  the  fuel  cell  works  at  the 
pre-established  optimal  operating  point.  For  these  reasons,  they  are 
vulnerable  to  faults  that  can  cause  the  stop  or  the  permanent  dam¬ 
age  of  the  fuel  cell.  To  guarantee  the  safe  operation  of  the  fuel  cell 
systems,  it  is  necessary  to  use  systematic  techniques,  like  the  recent 
methods  of  fault  tolerant  control  (FTC),  which  allow  to  increase 
the  fault  tolerance  of  this  technology  described  in  [1,2].  The  first 
task  to  achieve  active  tolerant  control  consists  of  the  inclusion  of  a 
fault  diagnosis  system  operating  in  real-time.  The  diagnosis  system 
should  not  only  allow  the  fault  detection  and  isolation  but  also  to 
the  fault  magnitude  estimation.  In  this  paper,  a  model-based  fault 
diagnosis  is  proposed  as  a  way  to  diagnose  faults  in  fuel  cell  systems. 
The  model-based  fault  diagnosis  is  based  on  comparing  on-line  the 
real  behaviour  of  the  monitored  systems  obtained  by  means  of  sen¬ 
sors  with  a  dynamic  model  of  the  same  simulated  system.  In  case  of 
a  significant  discrepancy  (residual)  is  detected  between  the  model 
and  the  measurements  obtained  by  the  sensors  the  existence  of  a 
fault  is  assumed.  If  a  set  of  measurements  is  available,  it  is  possible 
to  generate  a  set  of  residuals  (indicators)  that  present  a  different 
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sensitivity  to  the  set  of  possible  faults.  Analyzing  in  real-time  how 
the  faults  affect  to  the  residuals,  it  is  possible,  in  some  case,  to  iso¬ 
late  the  fault,  and  even  in  some  cases  it  is  also  possible  to  determine 
its  magnitude.  The  innovation  of  the  proposal  of  this  paper  is  based 
on  the  use  of  the  residual  fault  sensitivity  analysis  that  allows  to 
isolate  faults  that  otherwise  would  not  be  separable. 

The  structure  of  this  paper  is  the  following:  in  Section  2, 
the  foundations  of  the  proposed  fault  diagnosis  methodology  are 
recalled.  In  Section  3,  the  proposed  model-based  fault  diagnosis 
methodology  is  described.  In  Section  4,  the  PEM  fuel  cell  system 
used  to  illustrate  the  proposed  fault  diagnosis  methodology  is  pre¬ 
sented  with  the  fault  scenarios  that  can  appear.  Finally,  in  Section 
5,  the  application  results  of  the  proposed  methodology  of  diagnosis 
are  presented. 

2.  Foundations  of  the  fault  diagnosis  methodology 

2.1.  Model-based  fault  diagnosis 

The  methodology  of  fault  diagnosis  which  is  used  in  this  work  is 
mainly  based  on  classic  theory  of  model-based  diagnosis  described 
for  example  in  [3-6].  Model-based  diagnosis  can  be  divided  in 
two  subtasks:  fault  detection  and  fault  isolation.  The  principle  of 
model-based  fault  detection  is  to  check  the  consistency  of  observed 
behaviour  while  fault  isolation  tries  to  isolate  the  component  that 
is  in  fault. 

The  consistency  check  is  based  on  computing  residuals,  r(/<).  The 
residuals  are  obtained  from  measured  input  signals  u{k)  and  out- 
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puts  y{k)  using  the  sensors  installed  in  the  monitored  system  and 
the  analytical  relationship  which  are  obtained  by  system  modeling: 

r{k)  =  V{y{k),u(k))  (1) 


where  0  is  the  residuals  generator  function  that  depends  on  the 
type  of  detection  strategy  used  (parity  equation  [3]  or  observer  [7]). 
At  each  time  instance,  k,  the  residual  is  compared  with  a  threshold 
value  (zero  in  ideal  case  or  almost  zero  in  real  case).  The  threshold 
value  is  typically  determined  using  statistical  or  set-based  methods 
that  take  into  account  the  effect  of  noise  and  model  uncertainty 
[1  ].  When  a  residual  is  bigger  than  the  threshold,  it  is  determined 
that  there  is  a  fault  in  the  system;  otherwise,  it  is  considered  that 
the  system  is  working  properly.  In  practice,  because  of  input  and 
output  noise,  nuisance  inputs  and  modelling  errors  affecting  to  the 
considered  model,  robust  residuals  generators  must  be  used.  The 
robustness  of  a  fault  detection  system  means  that  it  must  be  only 
sensitive  to  faults,  even  in  the  presence  of  model-reality  differences 
[7]. 

Robustness  can  be  achieved  at  residual  generation  ( active )  or 
evaluation  phase  {passive).  Most  of  the  passive  robust  residual  eval¬ 
uation  methods  are  based  on  an  adaptive  threshold  changing  in 
time  according  to  the  plant  input  signal  and  taking  into  account 
model  uncertainty  either  in  the  time  or  frequency  domain  [8].  In 
this  paper,  a  passive  method  in  time  domain  has  been  proposed  for 
robust  fault  detection  in  time  domain,  where  the  detection  thresh¬ 
old  has  been  obtained  using  statistical  techniques. 

Robust  residual  evaluation  allows  obtaining  a  set  of  fault  signa¬ 
tures  (p{k)  =  [0i  (/<),  02  (/<), . . . ,  0n0(k)],  where  each  indicator  of  fault 
is  obtained  as  follows: 


0;(/<)  = 


r  o  if 
\i if 


|n(fc)|  <  ti 

n(k)  |  >  r,- 


(2) 


where  r ,•  is  the  threshold  associated  to  the  residual  rz(/<). 

Fault  isolation  consists  in  identifying  the  faults  affecting  the  sys¬ 
tem.  It  is  carried  out  on  the  basis  of  fault  signatures,  0,  (generated  by 
the  detection  module)  and  its  relation  with  all  the  considered  faults, 
/(/<)  =  {f(k)i,f2(k), . . .  ,fnf(k)}.  The  method  most  often  applied  is 
a  relation  defined  on  the  Cartesian  product  of  the  sets  of  faults 
FSM  c  0  xf  where  FSM  is  the  theoretical  signatures  matrix  [3].  One 
element  of  that  matrix  FSMy  will  be  equal  to  one,  if  a  fault  /•(/<), 
is  affected  by  the  residual  r2(/<).  In  this  case,  the  value  of  the  fault 
indicator  0z(/<)  must  be  equal  to  one  when  the  fault  appears  in  the 
monitored  system.  Otherwise,  the  element  FSMy  will  be  zero. 


which  is  a  transfer  function  that  describes  the  effect  on  the  residual, 
r,  of  a  given  fault/.  Sensitivity  provides  a  quantitative  information  of 
the  effect  of  the  fault  on  the  residual  and  a  qualitative  information 
in  their  sense  of  variation  (sign).  The  use  of  this  information  at  the 
stage  of  diagnosis  will  allow  separate  faults  that  even  presenting  the 
same  theoretical  binary  fault  signature,  presenting,  qualitatively  or 
quantitatively,  different  sensitivities. 

In  order  to  perform  diagnosis,  the  algorithm  will  use  the  the¬ 
oretical  signatures  FSMsensit.  FSMsensit,  as  any  FSM,  has  a  matrix 
structure  with  the  residual  sensitivity  in  the  row  and  the  faults  in 
columns,  each  value  of  this  matrix  will  be  notice  as  Sr.f..  Although 
sensitivity  depends  on  time  in  case  of  a  dynamic  system,  here  the 
steady-state  value  after  a  fault  occurrence  is  considered  as  it  is  also 
suggested  in  [3].  The  theoretical  value  of  Srjf.  describes  how  easily  a 
fault  will  cause  a  violation  of  the  threshold  of  the  ith  residual  since 
the  larger  its  partial  derivative  with  respect  to  the  fault /,  the  more 
sensitive  that  equation  is  to  deviations  of  the  assumption.  It  can 
computed  analytically  using  Eqs.  (1 )  and  (3)  or  in  simulation. 

In  order  to  perform  real  time  diagnosis,  the  observed  sensitivity 
S°fS  should  be  computed  using  the  current  value  of  the  residuals 

rvj 

rfk)  when  a  fault /(/<)  is  detected: 

m=s?p(k)  (4) 

3.  The  proposed  fault  diagnosis  methodology 

From  Eq.  (4),  it  can  be  seen  that  using  FSMsensit  in  real  time 
requires  the  knowledge  of  the  fault  magnitude  or  make  an  esti¬ 
mation  of  it.  To  solve  this  problem,  this  paper  attempts  to  design 
the  diagnosis  using  the  new  concept  of  relative  sensitivity  rather 
than  absolute  sensitivity  giving  in  Eq.  (3).  The  observed  relative  fault 
sensitivity  is  defined  as: 

crel.obs  =  = 

Sr,fj  n(k]fj(k)  r,(k)  1  ; 

which  corresponds  to  the  ratio  of  one  residue  rfk)  with  another  one, 
for  example  ri(/<).  Then,  the  relative  sensitivity  will  be  insensitive 
to  the  magnitude  of  that  unknown  fault. 

Using  the  concept  of  relative  sensitivity,  it  is  proposed  a  new 
FSM,  called  FSMsensit .rel,  which  corresponds  to  theoretical  fault  sig¬ 
nature  matrix  based  on  relative  sensitivities.  One  element  of  the 
theoretical  fault  signature  matrix  on  sensitivity  Srel’!.eo,  is  given  by: 

ri'uj 


2.2.  Fault  sensitivity  analysis 

The  isolation  approaches  presented  in  previous  section  uses  a 
set  of  binary  detection  tests  to  compose  the  observed  fault  signa¬ 
ture.  When  applying  this  methodology  to  dynamic  systems,  since 
they  may  exhibit  symptoms  with  different  dynamics,  the  use  of 
binary  codification  of  the  residual  produces  lose  of  information  [9]. 
This  can  be  the  origin  of  false  isolation  decisions,  especially  when 
some  detection  tests  have  a  transient  behaviour  (as  in  dynamic 
slow/delayed  systems)  in  response  to  the  faults.  Also,  in  complex 
systems,  some  faults  could  present  the  same  theoretical  binary  fault 
signature  not  allowing  fault  isolation.  In  both  cases,  it  is  possible  to 
use  other  additional  information  associated  with  the  relationship 
between  the  residuals  and  faults,  as  the  sign,  sensitivity,  order  or 
time  activation,  to  improve  the  isolation  results  [9]. 

In  this  work,  it  is  proposed  the  use  of  information  provided  by 
the  fault  residual  sensitivity  in  the  design  of  the  diagnosis  system 
in  order  to  increase  fault  isolability.  According  to  [3],  the  sensitivity 
of  the  residual  to  a  fault  is  given  by: 


crel.teo  S^'  foKjO/gfitjO 

Vifj  Sri/.  drm/dfjW  1  ; 

In  the  case  of  a  set  of  n  faults,  a  relative  fault  sensitivity  matrix 
FSMsensit. rel  should  be  used  as  the  one  shown  in  Table  1 . 

The  diagnostic  algorithm  will  be  used  to  assess  real-time 
observed  relative  sensitivities  using  Eq.  (5),  as  a  ratio  of  residuals, 
which  will  provide  a  vector  in  the  space  of  relative  sensitivities.  The 
vector  generated  will  be  compared  with  vectors  of  theoretical  faults 
stored  into  the  relative  sensitivity  matrix  FSMsensit. rel.  The  theo¬ 
retical  fault  signature  vector  with  a  minimum  distance  with  respect 


Table  1 

Theoretical  fault  signature  matrix  using  relative  sensitivity  with  respect  to  r i . 


/i 

fl 

fn 

r2/rl 

<-rel,teo 

Vi/i 

<-rel,teo 

r2r\h 

<-rel,teo 

r2rjn 

r3/n 

<-rel,teo 

r3rjt 

<-rel,teo 

r3rj2 

<-rel,teo 

r3rjn 

rmlr-i 

<-rel,teo 

rmrl/l 

<-rel,teo 

rmr\f2 

<-rel,teo 

rmrjn 

218 


T.  Escobet  et  al.  /  Journal  of  Power  Sources  192  (2009)  216-223 


to  the  fault  observed  vector  is  postulated  as  the  possible  fault: 

m\n{df](k), . . . ,  dfn(k)}  (7) 

where  the  distance  is  calculated  using  the  Euclidean  distance 
between  vectors: 

■W‘>  -  +<sS*’-0>2  <* 4 * * * 8> 

4.  Application  to  a  PEM  fuel  cell  system 

4.1.  Description  of  the  PEM  fuel  cell  system 

To  show  the  validity  of  the  proposed  model-based  fault  diagno¬ 
sis  approach  proposed  in  this  paper  when  applied  to  a  PEM  fuel  cell 
system  (PEMFC),  the  well  known  simulator  developed  by  Pukrush- 
pan  et  al.  [10]  is  used.  The  main  components  of  this  system  (Fig.  1) 
are  the  fuel  stack,  the  compressor,  the  inlet  and  outlet  air  man¬ 
ifold,  the  inlet  and  outlet  hydrogen  manifold  and  the  humidifier. 
The  air  supply  system  (compressor  and  air  collector)  has  as  pri¬ 
mary  objective  maintaining  a  constant  partial  pressure  of  oxygen 
in  the  cathode  outlet.  One  important  variable  is  the  oxygen  excess 
ratio  which  is  defined  as: 

y  _  ^C>2 , supplied  _  *02  ^air, supplied 
°2  =  Wq2,  reacted  =  JWo2n//c/4F 

where  x02  is  the  oxygen  molar  mass  fraction,  Wair  supplied  is  the  mass 
flow  rate  of  air  that  is  supplied  to  the  fuel  cell,  Mq2  is  the  oxygen 
molar  mass,  n  is  the  cell  number  in  the  stack,  IfC  is  the  stack  current, 
and  F  is  the  Faraday  number. 

The  aforementioned  model  is  a  control  oriented  model  for  con¬ 
trol  applications,  which  includes  the  transient  phenomenon  of  the 
compressor,  the  manifold  filling  dynamics  (both  anode  and  cath¬ 
ode),  reactant  partial  pressures,  and  membrane  humidity.  The  stack 
voltage  predicted  by  the  model  depends  on  load  current,  partial 


pressure  of  hydrogen  and  oxygen,  fuel  cell  temperature  and  the  con¬ 
tents  of  water  in  the  membrane.  Spatial  variations  are  not  included 
and  constant  properties  are  assumed  in  all  volumes.  Only  tempo¬ 
rary  variations  are  present.  The  model  also  assumes  that  the  inlet 
reactant  flows  in  the  cathode  and  in  the  anode  can  be  humidified, 
heated  and  cooled  instantaneously.  With  respect  to  the  considered 
dynamics,  the  model  neglects  the  fast  dynamics  of  the  electrochem¬ 
ical  reactions  (time  constant  of  10-19  s).  Temperature  is  treated  as  a 
constant  parameter  because  the  slow  behaviour  (time  constant  of 
102  s),  allowing  to  be  regulated  by  its  own  (slower)  controller.  The 
model  represents  a  75-kW  fuel  cell  system  with  381  cells. 

The  system  has  two  control  loops:  the  internal  loop  takes  the 
control  of  the  hydrogen  flow  and  the  external  loop  takes  the  con¬ 
trol  of  the  oxygen  excess  ratio  X0l ,  as  an  indirect  measure  to  control 
the  efficiency  of  the  PEMFC,  as  proposed  in  [10,11  ].  The  aim  of  the 
hydrogen  flow  control  is  to  minimize  the  differential  between  the 
anode  pressure  and  the  cathode  pressure.  The  regulation  of  is 
achieved  by  manipulating  the  compressor  motor  voltage  and  the 
output  cathode  valve.  The  stack  current  IfC ,  is  regarded  as  a  mea¬ 
sured  disturbance  to  the  system.  The  control  of  Xq2  is  achieved  using 
dynamic  matrix  control  (DMC),  a  control  technique  that  is  applied 
and  explained  in  detail  in  [11  ].  The  system  also  provides  measures 
of  the  compressor  current,  Icm,  and  its  speed,  &>cm.  Fig.  1  shows  an 
outline  of  the  PEMFC  system,  along  with  the  variables  available  for 
control  and  supervision  over  it. 

4.2.  Inclusion  of  the  faults  in  the  PEM  simulator 

In  this  paper,  the  PEM  fuel  cell  system  simulator  developed  in 
[10]  has  been  modified  in  order  to  include  a  set  of  typical  faults. 
The  faults  are  described  in  Table  2  and  the  description  of  how 
they  were  implemented  in  the  simulator  is  explained  in  the  fol¬ 
lowing. 

The  fault  f\  is  simulated  with  an  increment  A kv  in  the  com¬ 
pressor  constant  kv  and,  similarly,  the  fault  f2  is  simulated  with  an 
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Table  2 

Description  of  the  fault  scenarios. 

Fault  Description 

/i  Increase  in  the  friction  in  the  compressor  motor 

f2  The  compressor  motor  suffers  an  overheating 

/3  The  fluid  resistance  increases  due  to  water  blocking  the  channels  or 

flooding  in  the  diffusion  layer 

/4  Air  leak  in  the  air  supply  manifold 

/5  Increase  in  the  voltage  value  below  which  the  compressor  motor  does 

not  turn 

/6  Increase  in  the  stack  temperature  due  to  a  failure  in  the  temperature 

controller 


increment  A Rcm  in  the  compressor  motor  resistance  Rcm  •  Both  faults 
are  translated  in  a  change  in  the  compressor  torque  rcm  • 

?cm  =  p — ^  *  p — (Vcm  ~  (kv  +  A kv)C0cm)  (10) 

where  rjcm  is  the  motor  mechanical  efficiency,  kt  is  a  motor  constant, 
and  oocm  is  the  compressor  speed.  The  fault  f3  is  simulated  with  an 
increment  AWca,out  in  the  orifice  constant  of  the  cathode  output, 
kca, out,  which  produces  a  change  in  the  outlet  air  flow  in  the  cathode, 

Wca.out- 

Wca,  out  —  (kCa,  out  +  AkCa,out)(Pca  —  Prm)  (H) 

where  pca  is  the  cathode  pressure  and  pm  is  the  return  manifold 
pressure.  The  fault  f4  is  simulated  with  an  increment  A/<Sm,out  in 
the  supply  manifold  outlet  flow  constant  ksm, out,  which  is  translated 
into  a  change  in  the  outlet  air  flow  in  the  supply  manifold,  Wsm<0ut' 

Wsm  ,out  —  ( ksm  ,out  +  A/<sm,out)(P  sm  ~  Pca )  (12) 


Table  3 

Theoretical  fault  signature  matrix  FSM  using  binary  and  sign  information. 


Table  4 

Theoretical  fault  signature  matrix  FSMsensitsel. 


h 

h 

h 

k 

fs 

fs 

7*2 /*T 

1 

0.824 

0.118 

0.643 

0.036 

-0.221 

7*3 /^T 

0.854 

1 

0.197 

0.206 

0.039 

0.151 

7*4/7*! 

1 

0.937 

-0.128 

-0.134 

0.168 

-0.098 

expression: 

logio(Psat)  =  -1.69  X  10-10  T*+  3.85  x  1( )~7  T?c 

-3.39  x  10“4T|  +  0.14  Tfc-  20.92  (17) 

where  the  pressure  is  in  kPa  and  the  temperature  is  in  Kelvin.  The 
stack  temperature  also  affects  the  water  diffusion  coefficient  in  the 
membrane: 

D.-D),xp(24,6(5i5-5-3S¥))  (.8) 

where  is  a  constant,  which  depends  on  the  water  content  in  the 
membrane. 


where  psm  is  the  supply  manifold  pressure. 

The  fault  /s  is  simulated  with  an  increment  in  the  lower  volt¬ 
age,  Vcm  low,  that  the  controller  supplies  to  the  compressor  motor, 
a  boundary  that  also  influences  the  compressor  torque  in  Eq.  (10). 
The  fault  /6  is  simulated  with  an  increment  A TfC  in  the  stack  tem¬ 
perature  TfC ,  which  has  an  impact  on  the  open  circuit  voltage  of 
the  stack,  the  partial  of  gases,  the  relative  humidity,  and  the  water 
diffusion  coefficient  in  the  membrane.  The  open  circuit  voltage  E 
is: 


E  =  1.229  -  0.85  x  10“3(r/c  +  A Tfc  -  298.15) 

+  4.3085  x  10-5(r/c  + A7}c)[ln(pH2)  +  0.5  ln(po2)]  (13) 

where  pH2  and  p02  are  the  partial  pressure  of  hydrogen  and  oxy¬ 
gen,  respectively.  The  partial  of  gases  p*  in  the  anode  is: 

_  mi,anRi(Tfc  +  AT/c) 

Pi, an  —  i  / 

van 

where  the  subscript  i  is  either  H2 -hydrogen  or  v- vapor,  m  is  the 
molar  mass,  R  is  the  gas  constant,  and  Van  is  the  anode  volume.  The 
partial  pressure  of  gases  in  the  cathode  is: 


Pi,ca 


mi,caRi(Tfc  +  A7yc) 


(15) 


where  the  subscript  /  is  02 -oxygen,  N2 -nitrogen  or  v- vapor  and  Vca 
is  cathode  volume.  The  relative  humidity  </>j  is: 


<t>j  = 


Pv,j 

Psat(l/c  +  AT/-) 


(16) 


where  the  subscript  j  is  either  an- anode  or  ca-cathode  and  the 
saturation  pressure  of  vapor  is  calculatedusing  the  following 


4.3.  Residual  generation  and  fault  sensitivity  analysis 

The  set  of  measured  variables  are  Aq2,  03cm,  km  and  Vst.  Using 
these  variables  and  the  non-linear  model  presented  in  [10],  four 
residuals  can  be  derived: 

U  =  Ao2  —  ,  r2  =  Vst  ~  Vst,  r3  —  km  —  km, 

V4  =  03  cm  —  03  cm  (19) 

Using  the  PEMFC  fault  simulator  including  faults  described  in  Sec¬ 
tion  4.2,  it  has  been  determined  experimentally  if  the  faults  defined 
in  Table  2  affect  or  not  each  of  the  previous  residuals.  From  these 
results,  the  theoretical  binary  fault  signature  matrix  presented  in 
Table  3  can  be  derived.  It  can  be  noticed  that  all  the  considered 
faults  affect  all  the  residuals.  Thus,  the  faults  are  not  diagnosable. 
Even  taking  into  account  how  the  residual  sign  is  affected  by  faults, 
not  all  the  faults  are  diagnosable:  f\  cannot  be  distinguished  from 
f2  and/3  from/4. 

Alternatively  using  the  relative  fault  sensitivity  (5),  the  fault 
signature  table  FSMsensitsel  can  be  calculated.  The  values  of  this 
matrix  are  shown  in  Table  41 .  It  can  be  noticed  that  in  this  case  all 
the  considered  faults  are  isolable  since  the  following  condition  is 
satisfied: 

fora11  j+k  (20) 

This  can  be  seen  by  representing  the  values  of  the  theoretical  fault 
signatures  (see  Fig.  2)  in  the  three-dimensional  space  r2/ri,  r3/ri, 
and  r4/r! .  Since  there  is  no  overlapping,  the  six  faults  can  be  detected 
and  isolated  and  thus,  can  be  diagnosed. 


1  This  table  has  been  derived  from  PEMFC  in  the  normal  operating  point. 
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Fig.  2.  Theoretical  fault  signature  matrix  in  the  residual  space. 


4.4.  Implementation  of  the  fault  diagnosis  system  in  simulation 

The  implementation  of  the  fault  diagnosis  approach  proposed 
in  Section  3  is  done  in  the  MATLAB/SIMULINK  environment.  The 
fuel  cell  simulator  used  is  the  one  developed  in  [12]  at  the  Uni¬ 
versity  of  Michigan,  but  it  is  modified  to  allow  the  inclusion  of 
the  faults  described  in  Section  4.2.  In  this  simulator  (see  Fig.  3), 
it  is  possible  to  reproduce  any  of  the  faults  presented  in  Table  2. 
The  set  of  available  measurements  are  compared  with  their  pre¬ 
dicted  value  using  a  non-faulty  fuel  cell  model.  The  differences 


between  the  predicted  and  measured  values  generate  a  set  of  resid¬ 
uals  that  are  sent  to  the  fault  diagnosis  system.  When  a  fault 
appears,  the  residuals  that  are  sensitive  to  this  fault  take  a  value 
different  from  zero.  When  some  of  the  residual  values  cross  the 
detection  threshold,  the  fault  diagnosis  starts  reasoning  with  all 
the  violated  residuals.  The  reasoning  (described  in  Section  3)  is 
based  on  computing  the  minimum  distance  between  the  observed 
fault  signature  and  theoretical  one.  The  fault  that  approaches  the 
most  to  one  of  the  fault  signatures  is  the  one  indicated  as  a  possible 
fault. 

5.  Results 

In  order  to  evaluate  the  model-based  fault  diagnosis  methodol¬ 
ogy  proposed  in  this  paper,  the  fault  scenarios  and  fault  simulator 
presented  in  Section  4  have  been  used.  In  this  section,  results  of  two 
of  the  proposed  fault  scenarios  Cft  and  f2 )  are  presented.  The  results 
obtained  in  the  other  fault  scenarios,  although  they  have  not  been 
included  in  this  paper,  were  also  satisfactory  allowing  in  all  cases 
the  right  isolation  of  the  faults. 

5.1.  Fault  scenario  fo 

As  discussed  in  previous  sections,  the  fault  detection  is  based  on 
checking  at  every  time  the  difference  (residual)  between  the  sig¬ 
nals  monitored  by  a  sensor  and  its  estimation  using  the  detection 
model  (19).  Fig.  4(a)  shows  the  temporal  evolution  of  the  residuals 
and  the  detection  threshold  for  each  of  them.  The  fault  is  introduced 
into  the  system  at  time  50  s  creating  changes  in  its  internal  dynam¬ 
ics  and  sometimes  after  all  the  fault  signals  crosses  they  detection 
threshold  (dash  line)  activating  the  four  indicators  of  fault,  as  has 
been  defined  in  Eq.  (2).  Fig.  4(b)  illustrates  the  time  the  diagnosis 
system  takes  for  detecting  and  isolating  the  fault/!.  The  detection 


Fig.  3.  Implementation  of  the  fault  diagnosis  system  in  simulation. 
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Fig.  4.  (a)  Time  evolution  of  the  residual  corresponding  to  fault  fl,  (b)  Time  evolution  of  fault  detection  and  isolation  indicators  corresponding  to/i . 


subsystem  stores  the  fault  at  the  time  that  one  of  the  thresholds  is 
violated  by  any  of  the  residuals  and,  as  soon  as  it  is  detected  that 
a  fault  is  presented,  the  isolation  process  begins.  The  isolation  pro¬ 
cess  is  based  on  evaluating  the  distance  of  the  observed  relative 
fault  sensitivity  vector  to  the  theoretical  ones. 

Fig.  5(a)  shows  the  location  of  the  faults  in  the  space  of  ratios 
using  the  relative  fault  sensitivities  matrix,  described  in  Section 
3,  and  the  time  evolution  of  the  minimum  distance  between  the 
observed  and  theoretical  relative  fault  sensitivity  (7)  (drawn  in  con¬ 
tinuous  line). 

Fig.  5(b)  presents  the  Euclidean  distance  between  the  observer 
and  the  theoretical  sensitivity  fault  signatures  for  each  fault  (8).  It 
can  be  noticed  that  since  faulty  has  a  similar  fault  signature  as  f2 
(see  also  Table  4)  at  the  beginning  of  the  fault  isolation  process  f2 
is  the  fault  proposed  as  the  possible  fault  (since  presents  a  smaller 
distance  than/i ).  Flowever,  from  time  instant  82  s,/i  can  be  isolated. 


(a) 


It  is  seen  that  the  proposed  methodology,  after  some  isolation  time 
delay,  allows  isolate  the  true  fault. 

5.2.  Fault  scenario  f2 

Fig.  6(a)  shows  the  dynamic  evolution  of  residuals  when  the  fault 
f2  appears  in  the  system,  it  is  seen  that  residuals  have  the  same 
information  than /i.  Fig.  6(b)  illustrates  the  detection  and  isolation 
indicators  corresponding  to  fault /2.  The  process  of  isolation  before 
the  fault  f2  needs  more  processing  time  because  of  the  similarity 
between  the  fault  signatures  f\  and  f2. 

Fig.  7(a)  shows  the  location  of  faults  in  space  ratios  using 
the  matrix  of  theoretical  fault  signatures  and  the  time  evolution  of 
the  minimum  distance  corresponding  to  fault /2,  Fig.  7(b)  shows  the 
distance  of  the  observed  fault  signature  vector  compared  with  the 
theoretical  fault  signature  ones.  It  can  be  seen  that  the  fault  f2  can 


Fig.  5.  (a)  Time  evolution  of  the  minimum  distance  (7)  for  fault /i,  (b)  Time  evolution  of  the  distance  between  the  observed  and  theoretical  relative  fault  sensitivity  for  each 
fault. 
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Fig.  6.  (a)  Time  evolution  of  residual  corresponding  to  fault fj,  (b)  Time  evolution  of  fault  detection  and  isolation  indicators  corresponding  to/2. 
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Fig.  7.  (a)  Time  evolution  of  the  minimum  distance  for  fault/2,  (b)  Time  evolution  of  the  distance  between  the  observed  and  theoretical  relative  fault  sensitivity  for  each  fault. 


be  clearly  isolated,  because  always  corresponds  with  the  minimum 
distance  value. 

6.  Conclusions 

In  this  paper,  a  new  model-based  fault  diagnosis  methodol¬ 
ogy  based  on  the  relative  fault  sensitivity  has  been  presented  and 
tested.  An  advantage  of  this  new  methodology  is  that  it  does  not 
require  the  knowledge  of  the  fault  magnitude  to  provide  a  diagnos¬ 
tic.  Furthermore,  it  allows  isolate  faults  although  all  the  considered 
faults  affect  all  the  residuals  whenever  the  sensitivities  were 
different. 

To  prove  this  methodology,  a  PEM  fuel  cell  simulator  based  on 
the  model  presented  in  the  literature  has  been  developed.  The  sim¬ 
ulator  was  modified  to  include  a  set  of  possible  fault  scenarios 
proposed  in  this  work.  This  modified  simulator  allows  imposing 
a  determined  fault  scenario,  within  the  considered  set  of  faults  in 
the  PEMFC  and  analysing  its  behaviour.  All  the  simulated  faults 
have  been  tested  with  the  new  diagnosis  methodology,  which  has 
diagnosed  correctly  the  simulate  faults  in  contrast  with  other  well 
known  methodologies  using  binary  signature  matrix  of  analytical 
residuals  and  faults,  which  do  not  permit  to  isolate  the  complete 
set  of  faults. 
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